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ABSTRACT 

The material that formed the present-day Solar System originated in feeding zones in the inner 
Solar Nebula located at distances within ~ 20 AU from the Sun, known as the planet-forming zone. 
Meteoritic and cometary material contain abundant evidence for the presence of a rich and active 
chemistry in the planet-forming zone during the gas-rich phase of Solar System formation. It is a 
natural conjecture that analogs can be found amoung the zoo of protoplanetary disks around nearby 
young stars. The study of the chemistry and dynamics of planet formation requires: 1) tracers of 
dense gas at 100-1000 K and 2) imaging capabilities of such tracers with 5-100 (0.5-20 AU) milli-arcsec 
resolution, corresponding to the planet-forming zone at the distance of the closest star-forming regions. 
Recognizing that the rich infrared (2-200 jam) molecular spectrum recently discovered to be common in 
protoplanetary disks represents such a tracer, we present a new general raytracing code, RADLite, that 
is optimized for producing infrared line spectra and images from axisymmetric structures. RADLite 
can consistently deal with a wide range of velocity gradients, such as those typical for the inner regions 
of protoplanetary disks. The code is intended as a backend for chemical and excitation codes, and can 
rapidly produce spectra of thousands of lines for grids of models for comparison with observations. 
Such radiative transfer tools will be crucial for constraining both the structure and chemistry of 
planet-forming regions, including data from current infrared imaging spectrometers and extending 
to the Atacama Large Millimeter Array and the next generation of Extremely Large Telescopes, the 
James Webb Space Telescope and beyond. 

Subject headings: astrochemistry - radiative transfer - techniques: high angular resolution - planetary 
systems: protoplanetary disks 



1. INTRODUCTION 

The least processed components of the Solar System 
- comets, asteroids and Kuiper Belt Objects (KBOs) - 
are witnesses to an active Solar Nebula endowed with a 
rich chemistry. Much of the formation and early evolu- 
tion of the Solar System took place at small distances 
(perhaps within 20 AU) from the Sun where densities 
were high. Only at later stages did the solar system ex- 
pand to its present size (as defined by the Kuiper belt 
and Oort cloud), probably due to dynamical scattering 
processes involving planet esimals. It is thought that the 
KBOs were formed within 35 AU (Levison & Morbidelli 
2003), while the planets were all formed within 20 AU, 
Jupiter and Saturn within 10 AU (Tsiganis et al. 2005). 
Oort cloud comets are thought to have formed in the 5- 
10 AU region, and are composed of processed material 
that may have passed much closer to the Sun (Wooden 
et al. 2007). One group of important tracers of the Solar 
Nebula are primitive meteorites, including carbonaceous 
chondrites, since these can be brought into a laboratory 
setting and studied in great detail. Generally, such bod- 
ies formed in a relatively narrow range of radii from 2-4 
AU (Morbidelh et al. 2002). 

While the remnants of the formation of the Solar Sys- 
tem tell an exciting story, sometimes in great detail, they 
can never paint a full picture of planet formation in gen- 
eral because the evidence is ancient and hence often con- 
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taminated by processing and, perhaps more importantly, 
the Solar System is only one instance of a process that 
appears to be highly stochastic resulting in a wide range 
of planetary configurations (Udry & Santos 2007). For 
these reasons, the region spanning R = 0.1 — 20 AU 
in accretion disks around young stars can be considered 
the primary target for observational studies of exo-planet 
formation. The part of a protoplanetary disk within this 
range of radii is often putatively referred to as the planet- 
forming region. 

Given the angular sizes (5-100 milli-arcseconds) and 
surface temperatures (100-2000 K) of the planet-forming 
regions around T Tauri stars in the nearest star form- 
ing regions, they are most readily traced by molecular 
and atomic lines in the infrared - roughly 2-200 /im. 
Over the next decade, the mapping of the chemistry and 
structure of the planet-forming regions within hundreds 
of nearby protoplanetary disks, using infrared gas-phase 
signatures, can be expected to provide a critical comple- 
ment to the study of mature planetary systems around 
main-sequence stars. 

The first steps in this direction have already been 
taken, and an increasing number of infrared gas-phase 
tracers have been identified. The CO fundamental {v = 
1 — 0) ro- vibrational band at 4.7 /im is arguably the most 
accessible infrared molecular band. It has been used ex- 
tensively to trace gas in protoplanetary disks (Najita 
et al. 2003; Blake & Boogert 2004; Salyk et al. 2007), 
and CO fundamental line emission from disks has in a 
few cases been directly imaged using 8-10 class telescopes 
with adaptive optics (Goto et al. 2006; Pontoppidan et al. 
2008; van der Plas et al. 2009). Recently, Spitzer spec- 
troscopy has demonstrated that the 10-36 /im range of 




Fig. 1. — Sketch of the RADLite ray element subsamphng 
scheme. The ihustration shows the hne and continuum source Siy,a,s 
a function of radial velocity v{s), which in turn is a function of loca- 
tion s along the ray. Normally, the integral of the transfer equation 
is only evaluated at each end of the ray element, as indicated on 
the figure. However, lines with narrow local broadening may be 
completely missed by the integration. In such cases, the line is 
localized within the ray element and the integrand evaluated in a 
sufficient number of points across the line. 

disks is covered in molecular and atomic emission lines 
(Carr & Najita 2008; Salyk et al. 2008), dominated by 
H2O, OH, HCN, C2H2 and CO2, but also including im- 
portant atomic diagnostic lines such as the [Nell] transi- 
tion at 12.81 /im (Pascucci et al. 2007; Lahuis et al. 2007; 
Najita et al. 2009). In the visible, the [01] 6300 A line 
has also been used as an inner disk tracer of OH (Acke 
et al. 2005; Fedele et al. 2008; van der Plas et al. 2008). 

Among the central goals for kinematic imaging of 
planet-forming zones will be 1) a search for protoplanets 
and their dynamical interaction with their parent disks, 
2) a search for chemical processes inferred for the So- 
lar Nebula and 3) a general investigation of chemistry in 
planet-forming regions to better understand the initial 
conditions for planet formation. Key questions that can 
be answered include: If and how is water transported 
inside the snow line in protoplanetary disks, as was re- 
quired to form the Earth's oceans (see Encrenaz 2008, for 
a review)? Is it possible to draw direct connections be- 
tween protoplanetary disks and the Solar Nebula chem- 
istry inferred from primitive meteorites (Lyons & Young 
2005; Smith et al. 2009)? Are the inner holes and gaps 
observed in some protoplanetary disks associated with 
planet formation (Strom et al. 1989; Calvet et al. 2002; 
Lada et al. 2006; Brown et al. 2007)? Can protoplanets 
be detected dynamically through their interactions with 
the gaseous component of their parent protoplanetary 
disks? If so, what are their radial and mass distributions 
and how do these compare to the distributions of the 
population of mature exo-planets? Massive protoplanets 
are expected to form their own circumplanetary accre- 
tion disks (Canup & Ward 2002; Machida et al. 2008; 
Ayliffe & Bate 2009) - such as the disk that formed the 
Jovian moon system. Is it possible to detect circumplan- 
etary disks and use their kinematic signature to confirm 
the presence of a protoplanet? 

In this paper we present a new 2D ray-tracing code, 
called RADLite^, that is written specifically to rapidly 
and efficiently render images and spectra of lines in the 
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Fig. 2. — Unconvolved intensity image (I(x, y, v)) at Skms"-"^ for 
the 108 /im 22 1 li H2O line for the fiducial protoplanetary 
disk model viewed at an inclination of 45 degrees. For illustrative 
purposes, the maximal water abundance is lowered to 10~^ per 
hydrogen nucleus to keep the line optically thin. Due to the rela- 
tively long wavelength of this line, the disk is also has a dust optical 
depth of ~ 1, allowing the line image to show emission from both 
surfaces of the disk, resulting in line emission along two isovelocity 
contours. The white points indicate rays in which RADLite has 
detected line interactions. The far side of the disk is toward the 
top of the image. 
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Fig. 3. — Comparison of a ro- vibrational CO line rendered with 
(black curve) and without (grey histogram) the optimal gridding 
for the same computational time. This particular model with the 
optimal gridding turned off could be made to produce the same 
result as the optimal model, with a choice of a very fine global grid 
and a large number of rays, but at the cost of about a factor 10 in 
computational time. 

infrared, such as those tracing the planet-forming region, 
for a full gas+dust 2D density structure. Note, how- 
ever, that there is no fundamental reason why the code 
cannot also be applied to lines in the visible or (sub- 
) millimeter region. Most observational studies of infrared 
lines use very simple models, such as slab geometries or 
single layer thin disk calculations, to simulate lines from 
a growing database of high quality infrared line obser- 
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vations (e.g., Najita et al. 2003; Blake & Boogert 2004; 
Thi & Bik 2005; Brittain et al. 2007; Pontoppidan et al. 
2008). This is in sharp contrast to the relatively ad- 
vanced stage of radiative transfer modeling of rotational 
lines in the millimeter regime (Hogerheijde & van der 
Tak 2000; Pavlyuchenkov et al. 2007). With RADLite, 
we adopt a new philosophy: recognizing the very high 
complexity of inner disk chemistry and line excitation, 
we decouple these steps from the simulation and fitting 
of observational data. Because RADLite is fast, a user 
is able to calculate grids of models to determine best 
fitting absolute level population distributions in axisym- 
metric geometries. With these in hand, the application 
of, and comparison to, chemistry and excitation calcu- 
lations becomes a much simpler task. An example of 
such an application to the mid-infrared lines of water va- 
por from protoplanetary disks is presented in Meijerink 
et al. (2009). Here, the basic code is described in Section 
2 and the simulation of telescopic line imaging with the 
future generation of Extremely Large Telescopes (ELTs) 
is described in Section 3. 

As a demonstration, we apply RADLite to two exam- 
ples simulating ground-based high resolution (spectral 
and spatial) observations in Section 4: CO ro- vibrational 
emission from the fundamental {v = 1 — 0) band at 
4.7 /im (Section 4.2) and rotational water emission lines 
in the N-band at 10-14 /im (Section 4.3). These examples 
are used to demonstrate how RADLite simulates four 
common spectroscopic observables: line profiles, echel- 
lograms (also known as position-velocity diagrams or 2- 
dimensional slit images) , spectro-astrometry and integral 
field unit (IFU) image cubes. A 5th observable, interfer- 
ometric visibilities, is not explicitly simulated, but can 
easily be constructed from existing RADLite products, 
should the need arise. 

It is concluded, given that molecular lines in the near- 
to mid-infrared wavelength region are prime tracers of 
the dynamics and chemistry of the planet- forming re- 
gion, that the next generation of ground-based ELTs 
will be ideal for producing resolved line images. This 
requires an imaging spectrometer with resolving pow- 
ers of R = X/AX = 50,000 - 100,000. The Atacama 
Large Millimeter Array (ALMA), currently under con- 
struction, is also capable of high resolution molecular 
line imaging. This facility offers a tremendous improve- 
ment over previous millimeter arrays, and has sufficient 
spatial resolution to provide images with 1-2 AU resolu- 
tion of the planet-forming region in the continuum (Wolf 
& D'Angelo 2005). For line observations, ALMA will 
reach spatial resolutions of ^ 10 — 20 AU, or somewhat 
larger radii than those traced by infrared lines (Semenov 
et al. 2008). Further, some molecular species, such as the 
symmetric rotators CH4, CO2 and C2H2, available in the 
infrared cannot be reached with ALMA and vice versa. 
Perhaps the most notable example of a molecule gener- 
ally considered to be an infrared or space-based tracer 
is H2^0. Therefore, future observations of infrared lines 
tracing the inner disk (0.1-10 AU) will be highly comple- 
mentary to the (sub-) millimeter tracers to be imaged by 
ALMA, and powerful synergies between these facilities 
can be anticipated. 

2. RADIATIVE TRANSFER MODEL 
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Fig. 4. — Continuum-subtracted RADLite LTE renderings (solid 
curves) of the rotational ^^CO (3-2) and HCN (4-3) lines for a disk 
model designed to be as close as possible to that implemented for 
LkCa 15 by van Zadelhoff et al. (2001) (dashed curves). The differ- 
ence in strengths of the HCN lines can be explained by subthermal 
excitation in the outer disk in the van Zadelhoff et al. (2001) model, 
given the high critical density of this transition (1.2 x 10^ cm~^). 

2.1. A raytracer for infrared lines from protoplanetary 

disks 

RADLite is a general ray-tracing code that has been 
optimized for use with infrared lines in axisymmetric 
density distributions. It is intended to be used as a back- 
end to dust radiative transfer and line excitation calcula- 
tions. The continuum calculation should provide a dust 
temperature and the source function due to scattered 
photons in each grid point, as implemented in the con- 
tinuum Monte Carlo code RADMC (DuUemond & Do- 
minik 2004), while the line excitation calculation should 
provide level populations for a given molecular or atomic 
species. This construction is designed to take advan- 
tage of the growing availability of multi-dimensional line 
radiative transfer codes that need a tool to simulate ob- 
servational data (e.g., Woitke et al. 2009). In our imple- 
mentation, RADLite uses RADMC for the dust radiative 
transfer, which is recommended for use with RADLite. 
For the line excitation hetaSD (Poelman & Spaans 2005; 
Meijerink et al. 2008) is used in a ID-hlD configura- 
tion, but RADLite can interface with any axisymmet- 
ric line excitation code. RADLite does not assume that 
the grid spacing is regular in the radial or angular direc- 
tion. If the input chemical and/or excitation model has 
a fine gridding over a certain range of the angular coordi- 
nate, ^, perhaps to adequately sample steep gradients in 
abundance or temperature in the surface layer of a disk, 
RADLite can simply use that grid directly. If the user 
does not wish to carry out a detailed balance calculation 
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Fig. 5. — Spectrum of the fundamental ■'^^CO ro- vibrational band 
for a disk at inclination angles of 68 and 72 degrees, as rendered by 
RADLite, showing both emission and absorption components. The 
spectral resolving power is assumed to be very high (0.1 km s"-*^). 
This model illustrates that RADLite can efficiently, and within 
one instance, deal with spatial scales ranging from < 0.1 AU to 
> 100 AU as well as spectral scales ranging from < 0.1 km s"-*^ to 
the entire SED. 

for the level populations, RADLite can also render im- 
ages and spectra assuming that the level populations are 
those given by local thermodynamic equilibrium (LTE). 
In all cases, RADLite accepts molecular parameter in- 
put (levels energies, statistical weights and Einstein A 
values) directly from the HITRAN database (Rothman 
et al. 2005) or in the format of the Leiden Atomic and 
Molecular Database (LAMBDA, Schoier et al. 2005). In 
the following, we use parameters from HITRAN. 

RADLite is based on the "Long Characteristics" (LC) 
method of integration of the transfer equation in polar 
coordinates presented in Dullemond & Turolla (2000). 
LC integrations are performed by constructing a large 
number of parallel rays passing through a given density 
structure at different impact parameters. Each ray is 
then divided into a number of ray elements defined by 
the boundaries of the density structure grid elements. 
The formal transfer equation is integrated across each 
ray element, starting from infinity. For a dust contin- 
uum source function, this integration usually only re- 
quires that the density and temperature structures are 
well approximated by a linear function across the ray el- 
ement. RADLite includes the line emissivity function, 
jl, and opacity, a^, terms in the transfer equation 

^ = [at{R, 0) + al{R, 9, v)] 9, v) - /,] , (1) 

where the source function is: 

and V is the macroscopic velocity projected on the ray 
{R and 6 are the polar coordinates of the grid). 



The line source and opacity functions can vary on 
much smaller spatial scales than the continuum func- 
tions. Such cases typically result from a narrow intrinsic 
line profile coupled with large macroscopic velocity gra- 
dients. For a 1-dimensional model (radial motion only) 
with a monotonic velocity gradient, the problem could 
be solved by applying the Sobolov, or large velocity gra- 
dient (LVG), approximation. On a general axisymmetric 
grid, the LVG approximation cannot be assumed to hold, 
a priori, because a line may interact with itself at several 
different locations along a ray (an example of a situation 
where this happens is discussed in Section 4.2). For an 
LC integration that normally evaluates the integral of the 
transfer equation at the boundaries of the model grid, a 
small intrinsic line broadening relative to the global ve- 
locity gradient can result in such rapid changes of the 
source function that these are entirely contained within 
the end points of a ray element. This situation is illus- 
trated in Figure 1. The problem is essentially one of 
undersampling, and can be resolved if the disk model is 
defined on a very fine grid. However, such a brute force 
solution generally results in a great loss of computational 
efficiency because the vast majority of integration steps 
do not interact with a line. The result is very long pro- 
cessing times - more than 10 minutes per line on current 
workstations. Since an important goal of RADLite is to 
be able to render grids of models, with each model po- 
tentially containing thousands of lines (e.g. the forest 
of water lines in the infrared) , such inefficient raytracing 
would be unacceptable. 

How bad is this undersampling problem? To explain 
the viscous time scales of typical PP disks, the turbu- 
lent velocities are predicted to be small (Hartmann et al. 
1998). For low turbulence, intrinsic line widths are pre- 
dicted to be dominated by thermal broadening, which in 
some cold regions, such as the disk midplane, can be as 
low as tens of ms~^. In the inner disk, and for typical 
size grids used for the rendering of continuum images and 
spectra, the jump in velocities, Av, can be of order a few 
up to lOkms"^, typically at ~0.1 AU for a T Tauri star. 
In other words, the length, in velocity space, of a ray 
element may be 1-2 orders of magnitude larger than the 
width of the line. In these cases, the contribution from 
the line will simply not be included in the integration, 
causing an erroneous surface brightness along that ray. 
RADLite solves the problem by using an optimal form 
of grid refinement, while still allowing the global den- 
sity structure to be defined on a relatively coarse spatial 
grid - typically one designed for dust radiative trans- 
fer. For each integration step across a given ray element, 
the code checks for line crossings by assuming that the 
velocity projected along the ray is a smooth, linear func- 
tion. If a line crosses the ray within a ray element, a 
subdivision of the element is done in such a way that 
the line is integrated with an optimal number of samples 
(usually about 5 per FWHM of the line). The physical 
parameters, density, level population, etc., are calculated 
within a ray element using a simple linear interpolation 
between the values evaluated where the ray crosses the 
cell boundary. This approach is illustrated in Figure 1. 

2.2. Ray configuration 

The optimal sampling along a ray takes care of under- 
sampling in one dimension. However, proper rendering 



5 



of line spectra and image cubes also requires that the set 
of rays is arranged with a sufficiently fine spacing. For in- 
stance, in the case of a protoplanetary disk, infrared lines 
arise from a thin surface layer. This gives rise to the line 
emission being formed in a thin band along the isoveloc- 
ity curve in the intensity image at a given frequence or 
line velocity (see Figure 2). The isovelocity curve is a 
curve along which the projected velocity of the emitting 
disk surface is constant. A disk has two surfaces, and 
each surface will produce line emission along an isove- 
locity curve. If the disk is optically thick to continuum 
photons, only one isovelocity curve will be seen because 
dust blocks the observers view of the other surface. For 
illustrative purposes, the line in Figure 2 is chosen such 
that the disk has an optical depth close to unity, allowing 
the other surface to be seen. The smaller the intrinsic 
line broadening, the narrower the emission band. If the 
image is not sufficiently spatially resolved (i.e. if the pix- 
els are not small enough), the thin isovelocity band may 
easily be missed, again resulting in severe underestimates 
of line fluxes. The RADLite strategy to deal with this, 
without sacrificing efficiency, is to define a large number 
of closely spaced rays, but to only integrate those that 
actually interact with the line. All other rays are con- 
sidered "continuum rays" , and are integrated only once, 
assuming that the continuum emission does not change 
with frequency across the line. The rays are arranged in 
concentric circles with 150-600 rays per circle. The cir- 
cles are spaced according to the same radial logarithmic 
grid of the density structure. In the inner disk, an addi- 
tional 50-100 ray circles are added spanning inwards to 
the stellar surface in order to properly sample the inner 
rim. Line spectra are calculated by a simple surface in- 
tegral; the intensities of each ray circle are averaged, and 
the total intensity calculated by a 1-dimensional radial 
integral of the circle averages. The star is included as a 
special central ray. The user can specify whether to use 
a blackbody or a more complex stellar model. A typical 
ray configuration for optimal line rendering is shown in 
Figure 2. 

The user can modify the sampling of the rays, depend- 
ing on the application. For imaging, typically more rays 
are needed to avoid aliasing effects, i.e. artifacts caused 
by interpolating undersampled data from a circular grid 
to a rectangular grid. RADLite outputs the circular im- 
ages, i.e. the ray configuration and the associated inten- 
sities. The post-processing IDL script can then be used 
to resample the image on a rectangular mesh, with an 
arbitrary spatial sampling. The resampling uses TRI- 
GRID, a bilinear interpolation scheme in IDL. As with 
all interpolation schemes, undersampling can be a prob- 
lem, but the user can verify that flux is conserved by 
comparing to the surface integral of the original circular 
image. Typical interpolation errors are of order 1% in 
the results presented here. 

Including the enhancements discussed above, RADLite 
renders a spectrum and image- velocity cube of a line at 
Skms"-"^ resolution on a single 3 GHz Intel Xeon proces- 
sor core in 15-30 seconds. Line spectra can be rendered 
using ^ 50, 000 rays, while high quality images may re- 
quire up to 300,000 rays. 

2.3. Benchmark 
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Fig. 6. — Simulated astrometric spectra (SA(v), Equation 6) of 
the CO rovibrational ^' = 1 — P(8) line for the fiducial T Tauri star 
model. The data are simulated for a resolving power of Skms"-*^, 
as is offered by, e.g., CRIRES on the VLT. The top panel shows the 
fiux spectrum for a source at 125 pc, while the lower three panels 
show the astrometric spectra at three different slit position angles. 

There is currently no appropriate benchmark of 
RADLite, given the lack of other codes covering the 
same parameter space. However, it was confirmed that 
RADLite reproduces the result of RAYTRACE, the well- 
tested continuum raytracer associated with RADMC, to 
a fraction of a percent. This ensures that the dust ra- 
diative transfer is carried out correctly. In Fig. 3, we 
compare a line calculated with and without the RADLite 
subgridding and selective ray calculation. The ray grid- 
ding is adjusted such that the two lines take the same 



6 



number of CPU cycles to render. In practice, the dif- 
ference between results achieved with the RADLite sub- 
gridding, as compared to brute force calculations vary 
widely with line and underlying disk structure. Further- 
more, RADLite reproduces the results of the toy model 
fit to the SR 21 disk from Pontoppidan et al. (2008) (see 
Section 4.2). 

Finally, given that a significant volume of radiative 
transfer work is available for (sub) millimeter lines, a 
model was constructed to render the rotational lines of 
^^CO (3-2) and HCN (4-3) for a disk with the same pa- 
rameters as those used by van Zadelhoff et al. (2001) for 
the T Tauri star LkCa 15 based on a two-layer passive 
disk model (Chiang & Goldreich 1997). The main differ- 
ences is that our model is based on a full 2D dust model 
and assumes LTE level populations, while the van Zadel- 
hoff et al. (2001) does not. It can be noted that the crit- 
ical densities of the modeled lines are ~ 3.5 x 10^ cm~^ 
and ^ 1.2 x lO^cm"^, so it is more likely that HCN is 
out of LTE than CO (the midplane of the van Zadelhoff 
et al. (2001) disk has densities lower than the HCN (4- 
3) critical density beyond ~150 AU). This is reflected 
by a stronger RADLite LTE line. RADLite reproduces 
in detail the line strength and shape of the ^^CO (3-2) 
transition, which is known to be in LTE throughout the 
disk midplane. 

2.4. Caveats 

While RADLite is not limited to a particular geometry 
or velocity field, apart from the strict requirement that 
the density structure is limited to 2-dimensional polar co- 
ordinates, a number of caveats still apply. One property 
of the disk emission that is not modeled in a fully self con- 
sistent manner is scattering. Dust scattering of contin- 
uum photons is included, but is assumed to be isotropic. 
Energy emitted in the lines, however, is not able to scat- 
ter on dust. Including this would require a Monte Carlo 
approach, which would be extremely time consuming - 
especially if a large number of lines are to be simulated. 
Lines are also rendered assuming that there are no inter- 
actions with other lines. We find that this is generally 
a very good approximation for protoplanetary disks. In 
this case, line interactions will only happen when isove- 
locity contours overlap. For a Keplerian disk, this only 
occurs when the difference in line rest frequencies be- 
tween overlapping lines are of the order of the local line 
broadening or smaller, i.e. < lkms~^. In the general 
case, for different velocity fields and molecular species, it 
is possible to find cases where line overlaps are a prob- 
lem. One example are the Q branches of the mid-infrared 
HCN, CO2 and C2H2 bending rovibrational bands. The 
shift between transitions of different isotopologues (e.g., 
H2^0 to H2^0) is generally much too large to cause any 
problems with line overlaps. 

3. SIMULATING OBSERVATIONS 

The primary purpose of RADLite is to simulate obser- 
vations, and the package is therefore equipped with sim- 
ple post-processing scripts that can, given parameters for 
a virtual telescope, convert ideal images and spectra to 
three typical observational products: single-slit spectra, 
astrometric spectra, and image cubes (the last as would 
be observed with an idealized IFU). However, it is fore- 
seen that users typically will want to supply their own 



telescopic models to suit some required level of complex- 
ity of the telescope simulation. For this reason, this step 
is not considered an integral part of RADLite, but is pro- 
vided as an optional post-processing IDL script. In the 
following, the telescope simulation that is currently used 
is presented. Tools that are not yet developed include vis- 
ibility calculations for current and future infrared inter- 
ferometers, such as VLT-MATISSE (Lopez et al. 2008). 

3.1. Telescopic observables 

First, the way RADLite treats the basic observables; 
image cubes, echellograms, line profiles and spectro- 
astrometry is described. The simulations assume that 
observations are made with grating spectrometers (as is 
the case for almost all current mid-infrared astronomi- 
cal instrumentation). A grating spectrometer images a 
spectrally dispersed source on a 2-dimensional detector 
after the image of the source (a PSF for a point source) 
has been truncated by passage through a (approximately 
rectangular) slit with a length and a width, measured in 
angular units. This configuration results in a spectral im- 
age or echellogram that has a spatial dimension (along 
the slit) and a spectral dimension (across the slit). Some 
times more than one spectral order is imaged on the same 
detector in a so-called cross-dispersed configuration. It is 
also possible to place a number of slits next to each other 
and redirect the resulting spectra to different locations 
on the same detector. Such an instrument is known as 
a slicer-type IFU, and will produce a line image cube in 
which each pixel contains a spectrum (a so-called spaxel). 
Other types of IFUs exist but, for this paper, only a sheer 
is considered. A sheer IFU data cube can also be con- 
structed using a single-slit instrument by stepping the 
slit across the source, although this approach obviously 
results in a great loss in efficiency as compared to an 
IFU. 

The simulated observables are calculated as follows. 
The basic RADLite output is an intensity map as a func- 
tion of velocity offset, f, from a line center, /(x,7/,f), 
where the spatial coordinates, (x^y), are measured in 
angular units, assuming some distance to the source. A 
telescopic observation causes the intensity map to be con- 
volved with an image point spread function, PSF{x,y), 
and an instrumental line broadening profile, L{v), result- 
ing in an "telescopic" intensity map: 

/t(x, y, v) = [/(x, y, v) * PSF{x, y)] * L{v), (3) 

where * denotes convolution. 

The echellogram, as observed by a single slit grating 
spectrometer, is simply the telescopic intensity map in- 
tegrated over the entrance slit width. The slit is ori- 
ented at a position angle PA on the sky (negative ro- 
tation), with width, VF, such that the coordinates de- 
fined by the slit become, = xcos{PA) — ysm{PA) and 
rj = xsm{PA) + ycos{PA). In this system ^ and rj are 
coordinates across and along the slit, respectively. The 
echellogram becomes: 

nW/2 

E{rj,v)= / /T(r/,^,^)de (4) 

J-W/2 

The advantage of the echellogram is that it encompasses 
both spatial and spectral information. 
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Fig. 7. — Example of a VLT-CRIRES continuum-subtracted echellogram (E(r],v), Equation 4) of the P(8) CO line from the i; = 1 — 
fundamental band at 4.736 /xm toward SR 21. The line is marginally resolved in the spatial direction, with a blue-to-red peak offset of 
13 AU, as determined using spectro-astrometry (Pontoppidan et al. 2008). The right six panels show RADLite models of the line using the 
disk parameters of Pontoppidan et al. (2008), i.e. a 1.5 Mq, 6500 K, 3Rq star with a disk with an R = 6.5 AU gap, viewed at an inclination 
of 20 degrees. Both the observation and model have slits oriented at an angle of 10 degrees from the disk major axis. Otherwise, the disk 
is the fiducial model, with a CO abundance of 2 x 10~^ per hydrogen nucleus. The top panels is a model that assumes a gas-to-dust ratio 
of 1280 and is optically thick in the continuum at 4.7 /xm, the lower panels show a disk where the gas-to-dust ratio has been increased to 
2.6 X 10^, at which point the disk becomes optically thin to continuum photons at 4.7 //m (see the discussion in the text). The resolving 
power is 3kms~^ for both models and observation. The horizontal noise streaks is an artifact of the continuum subtraction due to the 
noise of the continuum frame being reproduced at each velocity. Hence, the streaks would also appear in real data. Note that the models 
assume diffraction limited imaging (FWHM = 122 milli-arcsec at 4.7)Lim), while the actual CRIRES observation probably has a significant 
seeing halo from incomplete AO correction. The point source sensitivity for the ELT observation is assumed to be 1.3 x 10~^ Jy rms for a 
1 hour observation, which is about 20 times better than that of VLT-CRIRES. 



The line profile is a higher level product, being the 
echellogram integrated over some chosen aperture along 
the slit {t]) with width, A. Often this integral is weighted 
by the signal strength in the r] direction, resulting in a so- 
called optimal extraction. However, since that approach 
can be difficult to quantify for extended sources, it is not 
used here. The line profile is: 



nA/2 

F{v)=Da X / E{r],v)d^ 

J -A/2 



(5) 



where Da is a constant factor, depending on the shape 
of the PSF, of order unity taking into account that A 
may not be large enough to contain all the source flux. 
If there is no noise added to the model, A can be chosen 
sufficiently large so that Da = 1- It is possible to also 
use the echellogram to calculate the astrometric (first 
moment) spectrum. Following Pontoppidan et al. (2008): 



SA{v) = Cax '-J'^ \ \ , 



(6) 



where Ca is a constant factor that approaches unity as 
A ^ oo, similar to Da- This spectrum is used to extract 
spatial information on scales smaller than the formal spa- 
tial resolution of the telescopic image. Note that if there 
is a strong continuum underlying the line, the astrometric 
spectrum of the line will be "diluted" by a factor depend- 
ing on the line-to-continuum ratio: (1 + i^continuum/^iine)- 



In other words, the measured astrometric spectrum as 
given in Equation 6 should be multiplied by this factor 
to determine the real, continuum-subtracted, extent of 
the line emission. 

3.2. PSF modeling 

In order to estimate the quality of telescopic line im- 
ages, an idealized diffraction PSF of an ELT is calculated 
and used to convolve model images. We use the analyt- 
ical approach of (Sabatke et al. 2005) relevant for pupil 
masks consisting of hexagonal segments. Note that we 
do not attempt to accurately model an actual telescope, 
since the various ELTs are still in the design phase. In 
particular, the influence on the pupil mask by the pres- 
ence of secondary mounts has been ignored, and the seg- 
ment alignment is assumed to be perfect. The adopted 
ELT mirror is inspired by that planned for the European 
ELT (E-ELT), and consists of 984 hexagonal segments, 
each 1.4 m wide, resulting in a 42 m aperture. 

3.3. Atmospheric transmission and correction for 
telluric lines 

A fundamental problem for ground-based infrared 
spectroscopy is the presence of strong absorption lines 
from molecules in the Earth's atmosphere (telluric lines). 
These lines are often the same transitions that may 
be targeted in astronomical objects, and protoplanetary 
disks in particular. If the telluric line is saturated, it 
is necessary to observe the target at a time of the year 
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Fig. 8. — Image cube (lT(x,y,v), Equation 3) of a simulated E-ELT IFU observation of the = 1 — P(8) CO line at 4.736 /im for 
SR 21 at an assumed distance of 125 pc, including an appropriate noise level, as discussed in the text. The far side of the disk is toward 
the top. The scale of the image is linear. This simulation can be compared to the simulated single slit observation in Figure 7. The first 
panel shows the continuum image obtained bluewards of the line. 90% of the continuum has been subtracted from the remaining panels to 
illustrate that the minimum required differential imaging contrast is only 1:10. 



when its radial velocity is large enough to shift the line 
away from its telluric counterpart. Some astronomical 
objects have large intrinsic radial velocities, but most 
nearby protoplanetary disks do not. Typically, shifts of 
order 20-30 km s~^ are required, which, by a coincidence 
of nature, are similar to the velocity of the Earth relative 
to the Sun. To simulate these circumstances, we use a 
transmission model calculated for Paranal in Chile, the 
location of the European Very Large Telescope, and as- 
sume an optimal, but generally achievable, velocity shift 
of 30 kms~^. 

Space-based infrared observatories, such as the James 
Webb Space Telescope (JWST), obviously do not re- 
quire the correction for telluric absorption. The Strato- 
spheric Observatory For Infrared Astronomy (SOFIA) 
has a significant atmospheric column above it for terres- 
trially abundant species such as CO2, CII4, or H2O, and 
hence still requires a correction for the telluric transmis- 
sion spectrum. 

3.4. Sensitivity 

An important component for predicting telescope per- 
formance is the sensitivity of the modeled observation. 
For a future ELT, predicting detailed sensitivities is 
clearly a futile task. For illustrative purposes, we use 
point source sensitivities of 125 /iJy (la in one hour) at 
M-band (4.6 /im) and 370 /iJy at N-band (10.5 /im) for re- 
solving powers of 100,000 and 50,000, respectively. These 
are values calculated using the online exposure time cal- 
culator for the E-ELT, available on the ESO website. The 
M-band sensitivity is only ~ 20 times better than that 
of VLT-CRIRES, which seems to be a conservative esti- 
mate. For instance, on the website for METIS, the high 



resolution mid-infrared spectrometer E-ELT concept, the 
point source sensitivities are reported to be ~10 times 
better at 4.6 /im. At N-band the gain seems to be signif- 
icantly higher compared to VLT-VISIR, according to all 
available estimates. Sensitivities for existing instrumen- 
tation are based on actual typical performance. 

4. MOLECULAR LINE IMAGING IN THE INFRARED 

A key property of ground-based telescopes with ex- 
tremely large apertures is their potential for delivering 
imaging with very high spatial resolution in the infrared 
- where the performance of adaptive optics systems is ex- 
pected to deliver nearly diffraction limited images with 
a resolution as high as 10-25 mas at 2-5 jam and 50-100 
mas at 10-20 /im. This corresponds to physical scales of 
2-10 AU in nearby star forming regions, which offer up 
to 1,000 potential targets including PP disks and other 
types of young stellar objects (YSOs) across the stel- 
lar mass range from the brown dwarf limit to young A 
stars (Evans et al. 2008). In this section, we demon- 
strate how RADLite can be applied to existing molecu- 
lar line spectra and explore what images of PP disks in 
various infrared molecular lines are expected to look like 
with a mid-infrared imaging spectrometer operating with 
spectral resolving powers of 50,000-100,000 mounted on 
an ELT. Several concepts for such ELT instruments 
are currently being developed. While the focus of the 
RADLite demonstration is on line imaging with ground- 
based facilities, there are other important facilities for 
mid-infrared spectroscopy, each covering unique aspects. 
The mid-infrared spectrometer on the JWST, MIRI, will 
provide highly sensitive medium resolution spectroscopy 
across the mid-infrared range (5-27 /im). For the bright- 



9 



est sources, the high resolution spectrometer on SOFIA, 
EXES, win enable spectrally resolved observations of a 
wide range of lines, not available from the ground, includ- 
ing the bending mode lines of water around 6 /im, in ad- 
dition to other important species such as CH4. RADLite 
is also intended to aid in the modeling of observations 
from these facilities. 

4.1. Fiducial disk model 

We base most of the examples on the fiducial disk 
model described in detail in a companion paper (model 
5 of Meijerink et al. 2009), but note the cases where 
departures from the fiducial model are made for illus- 
trative purposes, such as introducing an inner hole to 
simulate transition disks. Briefly, the central source is 
a M = 1 Mq , = 2 Rq star with an effective temper- 
ature of T = 4275 K, corresponding to a luminosity of 
1.4 L0 and an age of 2 Myr (Siess et al. 2000). The stel- 
lar input spectrum is a model atmosphere from Kurucz 
(1993). The disk has a gas mass of 10~^ M©, outer ra- 
dius of 120 AU, a Gaussian vertical density profile and 
a surface density that varies with radius as S ~ R~^. 
The inner rim is roughly located at the dust sublima- 
tion radius corresponding to 1600 K. The disk is assumed 
to have already experienced some grain growth, simu- 
lated by an enhanced gas-to-dust ratio of 12800, a low 
flaring index of H/R ~ and outer scale height 

^out/^out = 0.1. The total gas-to-solid ratio will of 
course still be ~100; the idea is to model the case where 
99% of the dust has already grown to centimeter or meter 
sizes, at which point they no longer contribute significi- 
cantly to the infrared opacity. Further, these large grains 
will likely have settled to the mid-plane, and will there- 
fore no longer contribute to the disk surface properties. 
The local line broadening is assumed to be dominated 
by the thermal component, corresponding to a low tur- 
bulent broadening, assumed to be a constant fraction of 
the local sound speed: ^turb = 0.03 Cg. The CO abun- 
dance is a constant 2 x 10~^ per hydrogen nucleus. The 
maximal water abundance (gas + ice) is taken to be high, 
3 X 10~^ per hydrogen nucleus (i.e. H+2H2) in shielded 
regions. As discussed in Meijerink et al. (2009), the in- 
nermost disk has the maximal water abundance, while 
the surface outwards of the mid-plane snow line, located 
at 1.7 AU, is strongly depleted in water down to a frac- 
tional abundance of 10~^^ per hydrogen nucleus in order 
to match observed mid-infrared water spectra in T Tauri 
disks. The level populations of the water vapor are calcu- 
lated in non-LTE using beta3D, while the CO is assumed 
to be in LTE. A complete, quantitative description of the 
RADLite input model, such as the density, temperature 
and abundance distributions, as well as a detailed treat- 
ment of the excitation of water vapor, can be found in 
Meijerink et al. (2009). 

4.2. The CO fundamental 

Given the large volume of observations of the CO fun- 
damental, it is expected that modeling of these lines will 
be a central application of RADLite. CO ro-vibrational 
lines are ubiquitous tracers of disk surfaces at ~ 1 AU 
(Najita et al. 2003; Blake & Boogert 2004; Brittain et al. 
2007; Salyk et al. 2007). Some disks are already being re- 
solved in CO ro-vibrational lines with AO-fed spectrom- 
eters on 8-10 m class telescopes (Goto et al. 2006). With 
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Fig. 9. — Top panel: A full range calculation of the infrared 
water spectrum from a typical protoplanetary disk (model 5 in 
Meijerink et al. (2009)). In addition, scaled (for clarity) model 
spectra convolved to the spectral resolution of various infrared ob- 
servatories are shown. For each observatory and instrument, rough 
sensitivities (~ lOcr in 1 hour) and spectral resolving powers are 
indicated. Bottom panel: A comparison of the model convolved 
to the Spitzer resolution (top curve) with Spitzer-IRS spectra of 
typical protoplanetary disks around solar-mass stars (from top to 
bottom, AA Tau, AS 205 and DR Tau, scaled for clarity, from Carr 
& Najita (2008); Salyk et al. (2008)). Note that the model only in- 
cludes water lines, not the Q branches from C2H2, HCN and CO2, 
which are prominent in the data. The observed spectra show lines 
of many more molecular and atomic species. 

the improvement in spatial resolution of an ELT by fac- 
tors of 3-5, the CO emission from many disks should be 
imaged over a correspondingly larger number of spaxels. 
Further, the CO fundamental can be used to demonstrate 
the ability of the code to simultaneously treat dynamic 
ranges of 4-5 orders of magnitude in both velocity and 
spatial dimension. Figure 5 shows a rendering of the 
^^CO band (lines from R(0)-P(14)) for the fiducial disk 
model viewed at high inclination (~ 70°). The models 
exhibit double-peaked emission line profiles typical for 
Keplerian disks. The lines are formed in the inner 1 AU 
of the disk. Superposed on the emission lines are very 
narrow absorption lines formed as continuum and line 
emission from the inner disk is absorbed by the outer 
disk - a configuration made possible by a fiaring disk ge- 
ometry. Such inclined disks have been used in the past 
to study disk material through absorption features (Pon- 
toppidan et al. 2005; Brittain et al. 2005; Lahuis et al. 
2006; Gibb et al. 2007). Two different models are pro- 
duced, for inclinations of 68 and 72 degrees, illustrating 
how rapidly the CO absorption changes with inclination; 
essentially no absorption is seen at inclinations below 68 
degrees, while the central continuum becomes very faint 
(<C O.lJy) at inclinations higher than 75 degrees. The 
prediction is that only ~5% of disks will have the appro- 
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priate inclination to show absorption from the outer disk. 
Quahtatively, the model spectrum at 72 degrees inclina- 
tion resembles that of HL Tau, a disk suggested to have 
a large gas-to-dust ratio in its upper layers, possibly to 
due to dust settling (Brittain et al. 2005). 

There are currently no high resolution IFU spectrom- 
eters that operate in the 3-20 /im region. However, one 
imaging technique that has been used for CO is spectro- 
astrometry (Pontoppidan et al. 2008; van der Plas et al. 
2009). Figure 6 shows astrometric spectra at three dif- 
ferent position angles (PAs) of the i = 72° P(8) CO ro- 
vibrational line also seen in Figure 5. At a PA=0°, the 
classical antisymmetric profile of a Keplerian disk ap- 
pears. The absorption from the outer disk shows up in 
the astrometry as a second antisymmetric profile - but 
with the opposite sign of that produced by the emission. 
This signature is due to the Keplerian rotation of the 
outer disk having a significantly varying radial compo- 
nent across the emission source (which consists of the 
central star, thermal emission from the innermost part 
of the disk, as well as scattering from the disk surface), 
and can therefore potentially be used to distinguish CO 
absorption from the disk, as opposed to CO absorption 
from an unrelated foreground cloud (which would have 
no astrometric signature). Also notable is a significant 
asymmetry seen in the off- axis angles 60 and 120°. For 
instance, at 60°, the blue-shifted lobe has a much smaller 
offset than the red-shifted lobe. This asymmetry is due 
to a combination of two effects. First, the structure of 
the inner rim in the disk causes an asymmetry of the 
continuum. When viewed at an inclined angle, the pho- 
tocenter of the continuum is shifted in the direction of 
the far side of the disk. Because spectro-astrometry is a 
relative measurement of a line compared to a continuum, 
this effect takes the appearance of a line asymmetry. Sec- 
ond, there is a true line asymmetry due to the fact that 
the disk is not perfectly fiat. Thus, the surface isovelocity 
contours are not symmetric around the disk major axis 
and the angle of a ray with respect to the disk surface 
normal is larger from the near side of the disk relative to 
the far side. The asymmetries are strong enough to be 
easily detectable with current spectro-astrometric facili- 
ties such as CRIRES on the Very Large Telescope ( VLT) , 
and can be effectively modeled with RADLite. 

The large apertures of the ELTs will be able to produce 
line images of planet- forming regions across many spax- 
els. The fundamental ro-vibrational band of CO will be 
the easiest imaging target. The expected improvement 
of a 42 m telescope relative to an 8.2 m telescope is il- 
lustrated in Figs. 7 and 8. The former shows an actual 
single-slit observation of a CO line from the protoplane- 
tary disk SR 21, which has an evolved planet-forming re- 
gion characterized by a large inner hole (a transitional or 
cold disk) (see Pontoppidan et al. 2008, for details). The 
observation is modeled with RADLite using the basic 
disk geometry derived by Pontoppidan et al. (2008), who 
determined that the CO gas is truncated within 6.5 AU. 
The additional applicability of RADLite is illustrated by 
simulating observations using two different gas-to-dust 
ratios to simulate the effect of grain growth. The low 
gas-to-dust ratio is 1280 while the high gas-to-dust ra- 
tio of 2.56 X 10^ has lowered the disk opacity sufficiently 
to make the disk optically thin to the continuum emis- 
sion at 4.7 /im. This produces two emission lobes in the 



position- velocity image, in contrast to the optically thick 
disk, which results in a single elongated structure. This 
difference is due to a limb brightening effect of the inner 
rim of the disk. For an optically thick inner rim, the 
surface intensity of the inner wall is constant at different 
ray impact parameters; the inner rim acts as an opaque 
"wall". An optically thin disk will result in a higher in- 
tensity at image locations where a ray is near tangential 
to the inner rim, i.e. where a ray samples the largest col- 
umn density of warm gas. While the observation hints at 
an optically thin disk, the greater sensitivity and higher 
spatial resolution of an ELT observation will enable such 
conclusions to be made with much greater confidence and 
for a much wider range of disks. 

It is likely that an ELT mid-infrared instrument will 
not be a single-slit facility, but an IFU. Figure 8 shows 
a RADLite simulation of the same disk as observed with 
a 42 m aperture telescope equipped with an IFU. In this 
case, the CO emission would be imaged across at least 50 
spaxels within 12 AU. It will thus be possible to search for 
kinematic structures in the disk directly related to planet 
formation, such as spiral waves and the radial density 
structure in the optically thin parts of the disk. 

4.3. H2O lines in the atmospheric N band window 

An important new window for high resolution spec- 
troscopy that will be opened by ELT class telescopes is 
the 7.5-13.9 /xm atmospheric N-band, as well as the 17- 
25 fim atmospheric Q-band. The only reason that these 
windows are currently not very productive for ground- 
based high resolution spectroscopic work, is that the sen- 
sitivity of instruments on 8-10 m class telescopes is not 
sufficient to reach the vast majority of protoplanetary 
disks, in particular disks around T Tauri stars now known 
to exhibit the richest molecular emission spectra. (See 
also recent work by Knez et al. (2009) and Najita et al. 
(2009) for discussions of the power and limitations of cur- 
rent high resolution N- and Q-band spectroscopy). The 
massive improvements in sensitivity offered by an ELT 
will put hundreds or thousands of protoplanetary disks 
within reach for high resolution imaging spectroscopy in 
the N- and Q-bands. 

There are many atomic and molecular species with 
strong transitions in the N- and Q-bands. One that 
has received considerable recent attention is the [Nell] 
12.8 /im line (Pascucci et al. 2007; van Boekel et al. 2009; 
Najita et al. 2009). Here, we show that a large number of 
rotational lines of water can be reached from the ground. 
These are lines that are known to be strong in many pro- 
toplanetary disks (Carr & Najita 2008; Salyk et al. 2008, 
and work in preparation by these authors). 

The strong lines and the associated high apparent 
molecular abundances in the disk surfaces within a few 
AU were not anticipated by the initial chemical mod- 
els - H2O abundances in the warm surface layers ap- 
pear to be of order nn^o/nn ~ 10~^ — 10~^, while 
many contemporary models (Glassgold et al. 2004; Gorti 
& Hollenbach 2008) predicted significantly less water - 
^H2o/^H ^ 10~^ — 10~^. Other modelers are now work- 
ing on explaining very high molecular abundances, in- 
cluding water (Glassgold et al. 2009; Woods & Willacy 
2009). 

Strong lines can be found across the entire 10-36 /im 
Spitzer high resolution spectral range and, by extension. 
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TABLE 1 

Strong N-band H2O lines suitable for ground- based observations 



Transition^ Wavelength [yLtm]^ Frequency [cm ^] Spitzer detected^ Strength^ Transmission^ 



V 


= 1- 





17o 17 




I81 18 


8.09583 


1235.20415 


N/A 


0.09 


clean 




V 


= 1- 





118 3 




129 4 


8.10842 


1233.28658 


N/A 


0.09 


clean 




V 


= 1- 





129 4 




13io 3 


8.19373 


1220.44594 


N/A 


0.09 


clean 




V 


= 1 - 





116 5 




127 6 


8.19441 


1220.34459 


N/A 


0.05 


clean 




p V 


= 1 - 





lie 6 




127 5 


8.20245 


1219.14720 


N/A 


0.05 


clean 




p V 


= 1 - 





I80 18 




19l 19 


8.23469 


1214.37420 


N/A 


0.08 


clean 




p V 


= 1- 





II5 7 




126 6 


8.24397 


1213.00716 


N/A 


0.07 


clean 




V 


= 1- 





128 5 




139 4 


8.26608 


1209.76296 


N/A 


0.06 


clean 




V 


= 1- 





127 6 




138 5 


8.32319 


1201.46258 


N/A 


0.05 


clean 




V 


= 0- 





136 7 




12l 12 


8.34329 


1198.56748 


N/A 


0.05 


clean 




V 


= 1- 





126 7 




137 6 


8.36579 


1195.34395 


N/A 


0.06 


clean 




V 


= 1 - 





138 5 




149 6 


8.42794 


1186.52920 


N/A 


0.05 


clean 




V 


= 1 - 





125 8 




136 7 


8.44520 


1184.10422 


N/A 


0.06 


clean 




V 


= 1- 





124 9 




135 8 


8.76006 


1141.54447 


N/A 


0.05 


clean 




V 


= 1- 





92 7 




IO5 6 


9.09356 


1099.67926 


N/A 


0.06 


clean 




V 


= 1- 





IO2 9 




II3 8 


9.16418 


1091.20494 


N/A 


0.06 


clean 




V 


= 0- 





146 9 




13i 2 


9.59209 


1042.52515 


N/A 


0.05 


clean 




V 


= 0- 





192 17 




I81 18 


9.88515 


1011.61860 


N/A 


0.09 


clean 




V 


= 0- 





14io 5 




137 6 


10.05642 


994.38945 


No 


0.05 


clean 




V 


= 0- 





177 10 




I64 13 


10.11319 


988.80730 


No 


0.05 


clean 




V 


= 0- 





178 9 




I65 12 


10.27482 


973.25273 


No 


0.05 


clean 




V 


= 0- 





13io 3 




127 6 


10.30317 


970.57472 


No 


0.05 


clean 




V 


= 0- 





I83 16 




17i 17 


10.42505 


959.22783 


No 


0.10 


clean 




V 


= - 





I69 8 




156 9 


10.47955 


954.23969 


No 


0.05 


clean 




V 


= 0- 





125 8 




llo 11 


10.54560 


948.26294 


No 


0.06 


clean 




V 


= 0- 





149 6 




136 7 


10.76434 


928.99316 


No 


0.07 


clean 




V 


= 0- 





127 6 




II2 9 


10.84439 


922.13542 


No 


0.08 


clean 




V 


= 0- 





1^6 9 




143 12 


10.85308 


921.39780 


No 


0.09 


clean 




V 


= 0- 





139 4 




126 7 


10.94112 


913.98290 


No 


0.08 


clean 




V 


= 0- 





194 15 




I83 16 


10.97414 


911.23324 


No 


0.07 


clean 




V 


= - 





I84 15 




17l 16 


10.98045 


910.70960 


No 


0.09 


clean 




V 


= - 





126 7 




111 10 


11.00168 


908.95248 


No 


0.10 


>30 kms- 


-1 


V 


= 0- 





158 7 




145 10 


11.02836 


906.75277 


No 


0.10 


>30 kms- 


-1 


p V 


= 0- 





173 15 




I60 16 


11.03376 


906.30900 


No 


0.05 


clean 




V 


= - 





172 15 
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11.03475 


906.22805 


No 


0.10 


clean 




p V 
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ll6 6 


11.23479 


890.09255 


No 
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clean 




V 


= 0- 
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11.25311 
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No 


0.10 


clean 




V 


= 0- 
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11.32409 


883.07302 


No 
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clean 




p V 
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II5 7 
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11.47773 
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No 


0.05 


clean 




V 


= 0- 





II9 2 
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11.55471 
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No 


0.10 


>30 kms- 


-1 


p V 


= 0- 
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11.56124 
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No 


0.05 


>30 kms- 


-1 


V 


= 0- 
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V 
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155 10 
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No 
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clean 




p V 


= 0- 
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clean 




p V 
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147 7 
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11.90026 


840.31789 
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clean 
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= 0- 





148 7 
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11.96812 
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No 
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clean 




p V 
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125 7 
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No 
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clean 
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175 13 
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12.24814 
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0.07 


clean 




V 
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II5 6 


12.27721 


814.51703 


Yes 


0.15 


>30 kms- 


-1 


V 


= 0- 





I64 13 




15i 14 


12.37566 


808.03788 


Yes 


0.15 


>30 kms- 


-1 


V 


= 0- 
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I63 14 


12.39625 
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Yes 


0.12 


clean 




p V 


= 0- 
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12.40705 
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Yes 


0.13 


clean 




V 


= 0- 
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IO5 6 


12.44483 


803.54638 


Yes 


0.17 


>30 kms- 


-1 


V 


= 0- 
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12.45346 
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Yes 


0.17 


>30 kms- 


-1 


p V 


= - 
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14o 14 


12.51462 


799.06551 


Yes 


0.15 


>30 kms- 


-1 


p V 


= - 





118 4 
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12.53831 
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Yes 


0.10 


clean 




V 


= - 
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12.59591 
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0.20 


>30 kms- 


-1 


p V 


= - 
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12.83197 
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No 


0.12 
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p V 


= - 
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12.89409 
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No 


0.15 
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-1 


p V 


= - 
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12.98575 
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No 


0.15 


>30 kms- 


-1 


V 


= - 





I65 12 
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13.03333 
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Yes 


0.12 


clean 




V 


= 0- 
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154 11 


13.13245 
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Yes 


0.09 


clean 




V 


= 0- 





153 12 




142 13 


13.29319 


752.26468 


Yes 


0.20 


>30 kms- 


-1 



^Bending mode vibrational quantum number. 
^Line parameters from HITRAN2004 (Rothman et al. 2005). 
^Subjective estimate limited by severe line blending. 
^Peak line flux/continuum flux. 

^Is a signiflcant Doppler shift required to observe the transition? 
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Fig. 10. — Simulations of ground-based observations of water lines from a typical protoplanetary disk. Top panel: model spectrum 
observed with a resolving power of R = A/ A A = 50, 000 with a conservative noise estimate for the E-ELT (3 times the theoretical estimate). 
Numerous strong lines (10-20% relative to the continuum in this particular model) are easily measured through the atmospheric transmission 
(second panel from the top). The black curve is the model spectrum with saturated (<50% transmission) parts of the telluric spectrum cut 
out. Third panel from the top: The same model spectrum, but now for the sensitivity (cr =100 mJy) of VLT-VISIR at A/AA = 30,000. 
Bottom panel: The observed Spitzer spectrum of AA Tau (lower black curve) relative to the same RADLite model (upper red curve) 
convolved to the resolution of the short-high module of Spitzer-IRS (A/AA = 600). All sensitivities assume roughly 1 hour of on-source 
integration time. Note that AA Tau is not a particularly bright target; with a continuum flux of 0.3 Jy at 10 /im, hundreds of protoplanetary 
disks will be brighter, some by an order of magnitude. 
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E-ELT 

HgO 108 3 -95 4 (12.87 ■ 



Inner edge at 0.07 AU 
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HgO 136 7 -123 10 (12.59 



Inner edge at 3 AU 




Inner edge at6AU 



-30 30 
Velocity [km/s] 

Fig. 11. — Simulated E-ELT continuum-subtracted echellograms 
of two N-band H2O lines for the fiducial non-LTE H2O model, both 
of which require velocity shifts of ~ SOkms"-*^ to shift the lines out 
of the telluric absorption line. The left line has an upper level 
energy of 3243 K, the right has an upper level energy of 3966 K. 
In the lower two rows of panels, models with a 3 and 6AU inner 
hole are also calculated. For clarity, the atmospheric transmission 
is not simulated. 

many strong lines are expected in the 6 /am water bend- 
ing mode as well as in the 40-200 /im range, in the regime 
of future space-borne observatories. Herschel-FACS will 
cover lines between 60 and 210 /im, but will not spec- 
trally resolve them. The first facilities with the potential 
to produce spatially and spectrally resolved data of rota- 
tional water lines from a significant sample of disks will 
be the ELTs. Figure 9 shows an overview of water lines 
along with near-term and future infrared facilities. The 
most sensitive ground-based window for observing the 
lines is the atmospheric A^-band window, which extends 
from 7.5 /im to 13.9 jam. Figure 10 shows a RADLite 
model spectrum of the fiducial model, based on observed 
Spitzer spectra^ of part of the N band. The model in- 
cludes atmospheric absorption for a dry site at 2600 m 
elevation. It also assumes that regions of the spectrum 
with less than 50% transmission cannot be calibrated, 
and these parts of the spectrum have been removed. 
As with most molecular infrared spectroscopy from the 
ground, it is necessary to observe the target at a time 
when the Doppler shift due to the relative motion of the 
Earth is large in order to move the target line away from 
its telluric counterpart. In this instance, a typical shift 
of 30kms~^ has been assumed. Finally, noise relevant 
for current mid-infrared high resolution spectrometers as 
well as for an instrument mounted on an ELT has been 
added to the spectra. It can be seen that for this particu- 
lar model, approximately 50 H2O lines will be observable 
from the ground under dry conditions (1-3 mm of precip- 
itable water vapor). Most line profiles can be measured 
in their entirety in one observing epoch, while a few of 
the lower excitation lines will require the combination 
of two epochs with opposite Doppler velocity shifts, i.e. 
separated by 6 months. The improvement in sensitivity 
of an ELT as compared to existing ground-based facil- 



ities clearly moves this type of observation from being 
very difficult to being straightforward for a large sam- 
ple of disks, given that the model example is of typical 
brightness for nearby T Tauri stars (0.35 Jy at 10 /im). 

Figure 11 shows simulations of II2O line imaging with 
the E-ELT, assuming a single-slit instrument. The fig- 
ure shows three different versions of the fiducial non- 
LTE model, the first with the inner disk truncated at 
the dust sublimation radius of 0.07 AU, the other two as- 
suming inner holes with radii of 3 and 6AU. The echel- 
lograms show that a 3 AU hole is easily spatially resolved. 
Two rotational lines are modeled: lOg 3 ^ 95 4 and 
136 7 1^3 10, with upper level excitation energies of 
3243 K and 3966 K, respectively (for comparison, the ex- 
citation energy of the CO ro- vibrational fundamental is 
~ 3100 K). The low excitation line can be detected out 
to 6 AU, while the higher excitation line drops below 
the detection limit between 3 and 6 AU. Astrometric 
spectra also detect the offsets for the 0.07 AU truncation 
radius. This shows that some lines are better tracers of 
certain radii than others due to a balance between exci- 
tation temperatures and geometric enhancement, i.e. the 
higher excitation line is the brightest at small radii but 
the faintest at large radii. 

In Table 1, an incomplete list of molecular lines most 
suitable for observations from the ground is given. The 
list is compiled based on line brightness (more than 5% 
relative to the continuum) and atmospheric transparency 
(more than 50% transmission over at least half the line, 
assuming 30kms~^ velocity shift relative to the atmo- 
sphere). Given that this selection is based on a model of 
one specific source, in practice, potentially thousands of 
molecular lines in 3-28 /im range will be detectable from 
the ground with an ELT equipped with a high resolution 
spectrometer. 

5. CONCLUSIONS 

We have presented a new axisymmetric ray tracer, 
called RADLite, for simulating telescopic line images, 
with specific emphasis on infrared lines formed in the 
planet-forming region of protoplanetary disks, and ap- 
plied it to mid-infrared CO and H2O emission. RADLite 
is optimized for rapidly rendering large numbers of lines 
for full 2D axisymmetric structures with arbitrarily large 
velocity gradients, including a detailed treatment of the 
dust radiative transfer. For instance, a full spectrum of 
~1000 water lines in the infrared (2-200 /im), can be ren- 
dered with a velocity resolution of 1 kms~^ in 1-2 hours 
on a single workstation. The code has very general appli- 
cations to chemical and excitation disk models as well as 
observations from infrared spectrometers on 8-10 m class 
ground-based telescopes, Spitzer-IRS^ Herschel-FACS^ 
future facilities such as the JWST, SOFIA, E-ELT, TMT 
and the Giant Magellan Telescope (GMT), in addition 
to (sub) millimeter facilities, such as ALMA. While not 
explicitly modeled in this paper, we also expect appli- 
cations to future infrared interferometer facilities, such 
as VLT-MATISSE. We find that a primary reason that 
infrared molecular spectroscopy of disks in the N-band 
(H2O and a wide range of other species) has not received 
much attention so far is due to a sensitivity deficit of 
roughly a factor of 10. This will be remedied by the 
ELT generation and the JWST, and we predict that 
this will become a highly active research area in the fu- 
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ture. The potential applications for RADLite reach be- 
yond protoplanetary disks. Any axisymmetric structure 
can be modeled. This could potentially include proto- 
planetary nebulae, the atmospheres of asymptotic giant 
branch stars, and active galactic nuclei. 
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